BABAR-PUB-06/032, SLAC-PUB-11902 Rev. Aug. 2006 
hep-ex/0606031 Phys. Rev. D74, 032005 (2006) 



Search for the decay B° K'^K'^K^ 

B. Aubert,^ R. Barate,^ M. Bona,^ D. Boutigny,^ F. Couderc,^ Y. Karyotakis/ J. P. Lees/ V. Poireau,"'^ 
V. Tisserand,! A. Zghiche/ E. Grauges,^ A. Palano,^ J. C. Chen ^ N. D. Qi,'^ G. Rong ^ P. Wang,"* Y. S. Zhu,^ 
G. Eigen,5 1. Ofte,^ B. Stugu,^ G. S. Abrams,^ M. Battaglia,*^ D. N. Brown,^ J. Button-Shafer,^ R. N. Gahn,^ 
E. Gharles,^ M. S. Gill,^ Y. Groysman,^ R. G. Jacobsen,^ J. A. Kadyk,^ L. T. Kerth,^ Yu. G. Kolomensky,^ 
G. Kukartsev,^ G. Lynch,^ L. M. Mir,*^ P. J. Oddone,^ T. J. Orimoto,'^ M. Pripstein," N. A. Roe,^ M. T. Ronan,*^ 
W. A. Wenzel,6 P. del Amo Sanchez/ M. Barrett/ K. E. Ford/ T. J. Harrison/ A. J. Hart/ C. M. Hawkes/ 
... S. E. Morgan/ A. T. Watson/ K. Goetzen/ T. Held/ H. Koch/ B. Lewandowski/ M. Pelizaeus/ K. Peters/ 
^ T. Schroeder/ M. Steinke/ J. T. Boyd/ J. P. Burke/ W. N. Cottingham/ D. Walker/ T. Cuhadar-Donszelmann/° 
S B. G. Fulsom/" C. Hearty/" N. S. Knecht/'' T. S. Mattison/" J. A. McKenna/'' A. Khan/^ P. Kyberd/^ 

^ . M. Salecm/i D. J. Sherwood/^ L. Teodorescu/^ V. E. Blinov/^ A. D. Bukin/^ V. P. Druzhinin/^ V. B. Golubev/^ 



< 



o 



X 



A. P. Onuchin/2 S. I. Serednyakov/^ Yu. I. Skovpen/^ E. P. Solodov/^ K. Yu Todyshev/^ D. S. Best/^ 
M. Bondioh/3 M. Bruinsma/^ M. Chao/^ S. Curry/^ I. Eschrich/^ D. Kirkby/^ A. J. Lankford/^ P. Lund/^ 
M. Mandelkern/3 R. K. Mommsen/^ W. Roethel/^ D. P. Stoker/^ S. Abachi/"* C. Buchanan/'' S. D. Foulkes/^ 
J. W. Gary/5 q Lo^g/^ B. C. Shen/^ K. Wang/^ L. Zhang/^ H. K. Hadavand/^ E. J. HiU/^ H. P. Paar/^ 
S. Rahatlou/s V. Sharma/^ J. W. BerryhiU/^ C. Campagnari/"^ A. Cunha/^ B. Dahmes/^ T. M. Hong/^ 
D. Kovalskyi/^ J. D. Richman/^ T. W. Beck/^ A. M. Eisner/^ C. J. Flacco/* C. A. Heusch/^ J. Kroseberg/* 
CN ; W. S. Lockman/8 G. Nesom/^ T. Schalk/^ B. A. Schumm/^ A. Seiden/^ P. Spradlin/^ D. G. WiUiams/^ 

' M. G. Wilson/8 J. Albert/^ E. Chen/^ A. Dvoretskii/^ D. G. Hitlin/^ I. Narsky/^ T. Piatenko/^ F. C. Porter/^ 
' A. Ryd/^ A. Samuel/9 R. Andreassen/o G. MancineUi/" B. T. Meadows/° M. D. SokolofF/" F. Blanc/^ 

O ■ P. G. Bloom/1 Chen/1 W. T. Ford/^ J. F. Hirschauer/i A. Kreisel/^ U. Nauenberg/i A. Olivas/i 

\0'. W. O. Ruddick/i J. G. Smith/1 ^ Ulmer/i S. R. Wagner/^ J. Zhang/i A. Ghen/^ E. A. Eckhart/^ 

O ■ A. Sofrer/2 W. H. Toki/^ R. J. Wilson/^ F. Winklmeier/^ Q. Zeng/^ D. D. Altenburg/^ E. Feltresi/^ 

A. Hauke/3 H. Jasper/^ A. Petzold/^ B. Spaan/^ T. Brandt/^ V. Klose/^ H. M. Lacker/'' W. F. Mader/"* 
R. Nogowski/4 J. Schubert/"* K. R. Schubert/'' R. Schwierz/" J. E. Sundermann/^ A. Volk/" D. Bernard/^ 

G. R. Bonncaud/'"^ P. Grenicr/^' * E. Latour/^ Gh. Thiebaux/^ M. Verderi/^ D. J. Bard/^ P. J. Glark/^ 
W. Gradl/6 F. Muheim/f^ S. Playfer/^ A. L Robertson/^ Y. Xie/^ M. Andrcotti/^ D. Bcttoni/"^ G. Bozzi/"^ 
R. Calabrese/"^ G. Gibinetto/'^ E. Luppi/"^ M. Negrini/'^ A. Petrella/"^ L. Piemontese/'^ E. Prencipe/'^ F. AnuUi/^ 
r-| R. Baldini-Ferroh/* A. Calcaterra/^ R. de Sangro/* G. Finocchiaro/^ S. Pacetti/^ P. Patteri/* L M. Peruzzi/*' t 
■ M. Piccolo/8 M. Rama/8 a. Zallo/^ A. Buzzo/^ R. Capra/^ R. Contri/^ M. Lo Vetere/^ M. M. Macri/^ 

M. R. Monge/9 S. Passaggio/^ C. Patrignani/^ E. Robutti/^ A. Santroni/^ S. Tosi/^ G. Brandenburg/" 
K. S. Ghaisanguanthum/" M. Morii/" J. Wu/" R. S. Dubitzky/i J. Marks/^ S. Schenk/i U. Uwer/i W. Bhimji/^ 
^ ! D. A. Bowerman/2 P. D. Dauncey/^ U. Egede/^ R. L. Flack/^ J. A. Nash/^ M. B. Nikohch/^ W. Panduro 

Vazquez/2 x. Ghai/^ M. J. Gharles/^ U. MaUik/^ N. T. Meyer/^ V. Ziegler/^ J. Gochran/'* H. B. Crawley/" 
L. Dong/" V. Eyges/" W. T. Meyer/" S. Prell/" E. L Rosenberg/" A. E. Rubin/" A. V. Gritsan/^ M. Fritsch/^ 
G. Schott/6 N. Arnaud/^ M. Davier/'^ G. Grosdidier/^ A. Hocker/"^ F. Le Diberder/"^ V. Lepeltier/^ 
A. M. Lutz/^ A. Oyanguren/^ S. Pruvot/"^ S. Rodier/^ P. Roudeau/"^ M. H. Schune/"^ A. Stocchi/"^ 
W. F. Wang/^ G. Wormser/^ C. H. Cheng/8 ^ Lange/8 D. M. Wright/8 C. A. Chavez/^ L J. Forster/^ 
J. R. Fry/9 g Gabathuler/^ R. Gamet/^ K. A. George/^ D. E. Hutchcroft/^ D. J. Payne/^ K. C. Schofield/^ 

C. Touramanis/9 A. J. Bevan/° F. Di Lodovico/*' W. Menges/° R. Sacco/° G. Gowan/i H. U. Flaecher/i 
D. A. Hopkins/1 P. S. Jackson/i T. R. McMahon/i S. Ricciardi/i F. Salvatore/i A. G. Wren/i D. N. Brown/^ 
C. L. Davis/2 Alhson/3 N. R. Barlow/^ R. J. Barlow/^ Y. M. Chia/^ C. L. Edgar/^ G. D. Lafferty/^ 
M. T. Naisbit/3 J. G. WiUiams/^ J. L Yi/^ C. Chen/" W. D. Hulsbergen/" A. Jawahery/" C. K. Lae/" 
D. A. Roberts/" G. Simi/" G. Blaylock/^ C. Dallapiccola/^ S. S. Hertzbach/^ X. Li/^ T. B. Moore/^ 
S. Saremi/5 H. Staengle/^ S. Y. Willocq/^ R. Cowan/^ G. ScioUa/^ S. J. Sekula/^ M. Spitznagel/^ F. Taylor/^ 
R. K. Yamamoto/6 H. Kim/^ P. M. Patel/^ S. H. Robertson/^ A. Lazzaro/8 V. Lombardo/8 F. Palombo/8 
J. M. Bauer/9 L. Gremaldi/® V. Eschenburg/" R. Godang/^ R. Kroeger/^ D. A. Sanders/^ D. J. Summers/^ 
H. W. Zhao/9 S. Brunet/o D. G6te/° P. Taras/° F. B. Viaud/" H. Nicholson/i N. Gavallo/^' i G. De Nardo/^ 

G. Sciacca/2 M. Baak/^ 



X 



2 



G. Ravcn,53 jj. L. Snock;'^-'' C. P. Jcssop,'^-* J. M. LoSecco,54 T. Allmcndinger,^^ q Bcnclli,-'^^ ^ Gan,^^ 
K. Honscheid,^^ D. Hufnagel,^^ P. D. Jackson,^^ H. Kagan,^^ R. Kass,^^ A. M. Rahimi,^^ R. Ter-Antonyan,^^ 

Q. K. Wong,55 N. L. Blount,^^ J. Brau,^^ R. Frey,^^ O. Igonkina,^^ M. Lu,^^ C. T. Potter,^^ R. Rahniat,^^ 
N. B. Sincv,56 q Strom,-'^^ j strubc,-'^^ g Torrcncc;"^^ F. Galcazzi,^^ A. Gaz," M. Margoni," M. Morandin," 
A. Pompili," M. Posocco," M. Rotondo,^^ F. Simonetto,^^ R. Stroili," C. Voci,^^ M. Benayoun,^^ J. Chauveau,^^ 
P. David,^^ L. Del Buono,^^ Ch. de la Vaissiere,^* O. Hamon,^^ B. L. Hartfiel,^^ M. J. J. John,^^ J. Malcles,^^ 
J. Ocariz,"^^ L. Roos,^^ G. Therin.^s P. K. Behera,^^ L. Gladney,^^ J. Panetta,^^ M. Biasini,™ R. Covarolli,60 
M. Pioppi,60 C. Angelini,6i G. Batignani,^! S. Bettarini,^! F. Bucci,<^i G. Calderini,^! M. Carpinelli,^! R. Cenci,^! 

F. Forti,^^ M. A. Giorgi,^^ A. Lusiani,^^ G. Marcliiori,^^ M. A. Mazur,^^ M. Morganti,^! N. Neri,^^ E. Paoloni,^^ 

G. Rizzo,^! J. Walsh,6i Haire,^^ ^ judd,*52 D. E. Wagoner,^^ J. Biesiada,^^ n. Daniclson,'^^ p Elmcr,63 

Y. P. Lau,63 G. Lu,^-'' J. 01scn,63 A. J. S. Smith,63 A. V. Telnov,63 F. Bellini,64 G. Gavoto,^^ A. D'Orazio,^^ D. del 
Re,*^^ E. Di Marco,'54 R. Faccini,^^ F. Ferrarotto,^'' F. Ferroni,^^ M. Gaspero,^"* L. Li Gioi,^'' M. A. Mazzoni,'''* 
S. Morganti,<54 G. Piredda,*^^ F. Polci,^^ F. Safai Telirani,^^ C. Voena,*^^ M. Ebert,*^^ H. Scliroder,^^ R. Waldi,<55 
T. Adyc,66 N. Dc Groot,^^ g Franck,^^ E. O. Olaiya,^^ F. F. Wilson,^^ R. Alcksan,^^ S. Emcry,^^ A. Gaidot,^'^ 
S. F. Ganzhur,^"^ G. Hamel de Monchenault,'^^ W. Kozanecki,*^^ M. Legendre,^'^ G. Vasseur,'^'^ Ch. Yeche,''^ 
M. Zito,<57 X. R. Clien,68 H. Liu,^^ W. Park,^^ M. V. Purohit,*^^ J. R. Wilson,^^ M. T. Allen,69 D. Aston,69 
R. Bartoldus,69 P. Bechtle,^^ N. Bergcr,^^ A. M. Boyarski,^'^ R. Claus,^^ J. P. Colcman,^^ M. R. Gonvcry,^^ 
M. Gristinziani,69 J. G. Dingfelder,'^^ J. Dorfan,^^ G. P. Dubois-Felsniann,^^ D. Dujmic,69 W. Dunwoodie,*^^ 

R. C. Field,69 T. Glanznian,69 S. J. Gowdy,^^ M. T. Graliani,69 V. Halyo,69 G. Hast,69 T. Hryn'ova,^^ 
W. R. Innes,69 M. H. Kelsey,^^ P. Kim,^^ D. W. G. S. Leith,^^ S. Li,^^ S. Luitz,^^ V. Luth,^^ H. L. Lynch,^^ 
D. B. MacFarlane,69 H. Marsiske,^^ R. Messner,^^ D. R. Muller,^^ G. P. O'Grady,'*^ y. E. Ozcan,69 A. Perazzo,^^ 
M. Perl,69 T. Pulliam,^^ B. N. Ratcliff,^^ a. Roodman,69 A. A. Salnikov,^^ R. H. Schindler,^^ J. Schwiening,^^ 
A. Snyder,69 J. Stclzcr,^^ D. Su,^^ M. K. Sullivan.^s K. Suzuki,^^ S. K. Swain,'^'' J. M. Tliompson,69 J. Va'vra,'^^ 

N. van Bakel,^'^ M. Weavei,^^ A. J. R. Weinstein,^^ W. J. Wisniewski,^^ M. Wittgen,^^ D. H. Wright,'^^ 
A. K. Yarritu,69 K. Yi,''^ C. C. Young,^^ P. R. Burchat,™ A. J. Edwards,™ S. A. Majewski,™ B. A. Petersen,™ 

C. Roat,™ L. Wilden,™ S. Alimed,^! M. S. Alam,^! R. Bula,^i J. A. Ernst,^! V. Jain,^! B. Pan,^! 
M. A. Saccd,'^! F. R. Wapplcr,^i S. B. Zain,'^i W. Bugg,^^ ]vj Krishnamurthy,'^^ g jyj Spanicr,^^ ^ Eckmann,'^^ 
J. L. Ritchie,™ A. Satpathy,™ C. J. Schilling,™ R. F. Schwitters,™ J. M. Izen,^" i_ Kitayama,^'' X. C. Lou,''^ 
S. YeJ^ F. Bianchi,™ F. Gallo,™ D. Gamba,™ M. Bomben,™ L. Bosisio,™ C. Cartaro,™ F. Cossutti,™ 

G. Delia Ricca,™ S. Dittongo,™ S. Grancagnolo,™ L. Lanccri,™ L. Vitalc,™ V. Azzolini,'^'^ F. Martinez- Vidal,'^^ 
Sw. Banerjee,™ B. Bhuyan,™ C. M. Brown,™ D. Fortin,™ K. Hamano,™ R. Kowalewski,™ I. M. Nngont,™ 

J. M. Roney,™ R. J. Sobie,™ J. J. Back,™ P. F. Harrison,™ T. E. Latham,™ G. B. Mohanty,™ M. Pappagallo,™ 

H. R. Band,80 X. Ghen,^" B. Chcng,^^ S. Dasu,8° M. Datta,8° K. T. Flood,«o J. J. Hollar,^^ P. E. Kutter,^^ 
B. Mellado,80 A. Mihalyi,80 Y. Pan,«o M. Pierini,80 R. Prcpost,^" S. L. Wu,^° Z. Yu,«° and H. Neal^i 

(The BABAR GoUaboration) 

^ Laboratoire de Physique des Particules, F-74941 Annecy-le-Vieux, Prance 

^Universitat de Barcelona, Facultat de Fisica Dept. ECM, E-08028 Barcelona, Spain 
^ Universita di Ban, Dipartimento di Fisica and INFN, 1-70126 Bari, Italy 
^Institute of High Energy Physics, Beijing 100039, China 
''' University of Bergen, Institute of Physics, N-5007 Bergen, Norway 
^Lawrence Berkeley National Laboratory and University of California, Berkeley, California 94720, USA 
^University of Birmingham, Birmingham, B15 2TT, United Kingdom 
^Ruhr Universitat Bochum, Institut fiir Experimentalphysik 1, D-4i780 Bochum, Germany 
^University of Bristol, Bristol BS8 ITL, United Kingdom 
University of British Columbia, Vancouver, British Columbia, Canada V6T IZl 
'^'^ Brunei University, Uxbridge, Middlesex UB8 3PII, United Kingdom 
^^Budker Institute of Nuclear Physics, Novosibirsk 630090, Russia 
University of California at Irvine, Irvine, California 92697, USA 
University of California at Los Angeles, Los Angeles, California 90024, USA 

University of California at Riverside, Riverside, California 92521, USA 
^''University of California at San Diego, La Jolla, California 92093, USA 
''^University of California at Santa Barbara, Santa Barbara, California 93106, USA 
'^University of California at Santa Cruz, Institute for Particle Physics, Santa Cruz, California 95064, USA 
'^California Institute of Technology, Pasadena, California 91125, USA 
^'^ University of Cincinnati, Cincinnati, Ohio 45221, USA 
^' University of Colorado, Boulder, Colorado 80309, USA 
Colorado State University, Fort Collins, Colorado 80523, USA 



3 



Universitat Dortmund, Institut Jiir Physik, Dortmund, Germany 

Technische Universitat Dresden, Institut fur Kern- und Teilchenphystk, D-01062 Dresden, Germany 
^^Ecole Poly technique, LLR, F-91128 Palaiseau, France 
University of Edinburgh, Edinburgh EH9 3JZ, United Kingdom 
Universita di Ferrara, Dipartimento di Fisica and INFN, I-44IOO Ferrara, Italy 
^^Laboratori Nazionali di Frascati dell'INFN, I-OOO44 Frascati, Italy 
^^Universita di Genova, Dipartimento di Fisica and INFN, I-I6I46 Genova, Italy 
'^"Harvard University, Cambridge, Massachusetts 02138, USA 
Universitat Heidelberg, Physikalisches Institut, Philosophenweg 12, D-69120 Heidelberg, Germany 
Imperial College London, London, SW7 2AZ, United Kingdom 
University of Iowa, Iowa City, Iowa 52242, USA 
^■^lowa State University, Ames, Iowa 50011-3160, USA 
Johns Hopkins University, Baltimore, Maryland 21218, USA 
^^Universitat Karlsruhe, Institut fiir Experimentelle Kernphysik, D-76021 Karlsruhe, Germany 
Laboratoire de I'Accelerateur Lineaire, IN2P3-CNRS et Urnversite Paris-Sud 11, 
Centre Scientifique d'Orsay, B.P. 34, F-91898 ORSAY Cedex, France 
Lawrence Livermore National Laboratory, Livermore, California 94550, USA 
University of Liverpool, Liverpool L69 7ZE, United Kingdom 
Queen Mary, University of London, El 4NS, United Kingdom 
University of London, Royal Holloway and Bedford New College, Egham, Surrey TW20 OEX, United Kingdom 
University of Louisville, Louisville, Kentucky 40292, USA 
'^^ University of Manchester, Manchester M13 9PL, United Kingdom 

University of Maryland, College Park, Maryland 20742, USA 
'^^ University of Massachusetts, Amherst, Massachusetts 01003, USA 
'^^Massachusetts Institute of Technology, Laboratory for Nuclear Science, Cambridge, Massachusetts 02139, USA 
^''McGill University, Montreal, Quebec, Canada H3A 2T8 
^^Universita di Milano, Dipartimento di Fisica and INFN, 1-20133 Milano, Italy 
'^^ University of Mississippi, University, Mississippi 38677, USA 
^°Universite de Montreal, Physique des Particules, Montreal, Quebec, Canada H3C 3J7 
Mount Holyoke College, South Hadley, Massachusetts 01075, USA 
^^Universita di Napoli Federico II, Dipartimento di Scienze Fisiche and INFN, 1-80126, Napoli, Italy 
NIKHEF, National Institute for Nuclear Physics and High Energy Physics, NL-1009 DB Amsterdam, The Netherlands 

University of Notre Dame, Notre Dame, Indiana 46556, USA 
^^Ohio State University, Columbus, Ohio 43210, USA 
University of Oregon, Eugene, Oregon 97403, USA 
Universita di Padova, Dipartimento di Fisica and INFN, 1-35131 Padova, Italy 
Universites Paris VI et VII, Laboratoire de Physique Nucleaire et de Hautes Energies, F-75252 Paris, France 
University of Pennsylvania, Philadelphia, Pennsylvania 19104, USA 
''^ Universita di Perugia, Dipartimento di Fisica and INFN, 1-06100 Perugia, Italy 
^' Universitd di Pisa, Dipartimento di Fisica, Scuola Normale Superiore and INFN, 1-56127 Pisa, Italy 
Prairie View A&M University, Prairie View, Texas 7744^, USA 
Princeton University, Princeton, New Jersey 08544, USA 
^'^ Universita di Roma La Sapienza, Dipartimento di Fisica and INFN, 1-00185 Roma, Italy 
''^'Universitat Rostock, D-18051 Rostock, Germany 
Rutherford Appleton Laboratory, Chilton, Didcot, Oxon, 0X11 OQX, United Kingdom 
^'^DSM/Dapnia, CEA/Saclay, F-91191 Gif-sur-Yvette, France 
University of South Carolina, Columbia, South Carolina 29208, USA 
^^Stanford Linear Accelerator Center, Stanford, California 94309, USA 
^Stanford University, Stanford, California 94305-4060, USA 
'^^ State University of New York, Albany, New York 12222, USA 
''^ University of Tennessee, Knoxville, Tennessee 37996, USA 
''^University of Texas at Austin, Austin, Texas 78712, USA 
University of Texas at Dallas, Richardson, Texas 75083, USA 
''^ Universita di Torino, Dipartimento di Fisica Sperimentale and INFN, 1-10125 Torino, Italy 
''^ Universitd, di Trieste, Dipartimento di Fisica and INFN, 1-34127 Trieste, Italy 
'^''IFIC, Universitat de Valencia- CSIC, E-46071 Valencia, Spain 
''^^ University of Victoria, Victoria, British Columbia, Canada V8W 3P6 
Department of Physics, University of Warwick, Coventry CV4 7AL, United Kingdom 
^"University of Wisconsin, Madison, Wisconsin 53706, USA 
^' Yale University, New Haven, Connecticut 06511, USA 

We present the first search for the decay — > K"K"K" using a data sample of 232 million BB 
pairs. We find no statistically significant evidence for the non-resonant component of this decay. 
Our central value for the branching fraction, assuming the true Dalitz distribution is uniform and 
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excluding the 4> resonance, is B{B° K°sK°sK°) = (2.4t|^ ± 0.6) x lO"*^ where the errors arc 
statistical and systematic, respectively. We set a single-sided Bayesian upper limit of B{B^ —> 
KgKgKl) < 7.4 X 10~^ at 90% confidence level using a uniform prior probability for physical 
values. Assuming the worst-case true Dalitz distribution, where the signal is entirely in the region 
of lowest efficiency, the 90% confidence level upper limit is B{B° -> K'^K^KI) < 16 x 10"^. 

PACS numbers; 13.25.Hw, 12.15.Hh, ll.30.Er 



I. INTRODUCTION 

Measurements from the BABAR and Belle experiments 
have confirmed the Cabibbo-Kobayashi-Maskawa quark 
mixing matrix (CKM) mechanism [1] as the dominant 
source of CP violation in flavor-changing weak interac- 
tions [2]. In the Standard Model, the mixing- induced, 
time-dcpcndcnt CP asymmetry in decays from b 
sqq penguin transitions should be the same as the pre- 
cisely measured CP asymmetry in charmonium-ii'° de- 
cays, namely sin 2/3, to within a few percent. New physics 
from higher mass scales may contribute to the loop in the 
penguin diagram, which could significantly alter the CP 
asymmetry in penguin-dominated B decays [3]. Initial 
CP asymmetry measurements in 6 — ^- sqq penguin B de- 
cays have suggested a possible violation of this test of the 
Standard Model [4], [5], [6]. 

The b sqq penguin decays fall into two categories. 
If the qq can be uu, a CKM-suppressed tree-level b 
u transition can contribute to the decay in addition to 
the dominant b sqq penguin. This introduces some 
uncertainty in the Standard Model prediction of the CP 
asymmetry, since the b ^ u and the penguin amplitudes 
have different weak phases. On the other hand, decays 
such as B" K",K",K", [5] and B" ^ ^s^s^" arc 
purely b sss penguin transitions and they can only 
include b ^ u decay amplitudes through rescattering, 
thus the Standard Model uncertainty on the predicted 
CP asymmetry for these decays is generally smaller [7] . 

It has been noted that three-body B^ decays of the 
form B° PPP', where P and P' arc spin-0 CP 
eigenstate neutral particles, are CP eigenstates[8], thus 
the CP asymmetry in B° K'^K^K'l is not diluted, 
as is generally the case in three-body i?" decays. The 
B" KgK'^Kl would be a valuable addition to under- 
standing the b sqq penguin CP asymmetry anomaly, if 
the branching fraction is large enough. The resonant (pK'^ 
contribution to this final state, neglecting interference 
effects, would not yield a sample of signal events large 
enough to make an interesting CP asymmetry measure- 
ment. In the analysis described below, the expected and 
observed number of K'^K^Kl events with a K^K^ invari- 
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ant mass in the region of the (j) resonance is about eight 
events. Prior to our search, the non-resonant component 
of the K'^K'^K'^ final state had not been experimentally 
investigated. This search was motivated by the possi- 
bility of discovering a large non-resonant contribution to 
the K'^K'^KI final state. Large non-resonant S-wave am- 
plitudes have been found in the Dalitz amplitude anal- 
ysis of B+ K+K-K+ [9] and ^ K+ K' [10]. 
The branching fraction was recently estimated to be [7] 
B{B^ - KtKtKl) = {^.2?,t\f,tr^^V£) X 10-« in- 
eluding the (j) resonance. The difference between the 
upper and lower limits of this estimate is substantial. 
Another prediction, based on isospin and Bose symme- 
try [11], gives a somewhat small branching fraction of 
B{B^ ^ K^,K^,Kl) = (1/3) B{B^ - K^sK^K^s) « 
2 X 10-*^. 



II. THE BABAR DETECTOR AND DATASET 

The resiilts presented here arc based on data col- 
lected with the BABAR detector [12] at the PEP-11 asym- 
metric e+e~ collider [13] located at the Stanford Lin- 
ear Accelerator Center. An integrated luminosity of 
211 fb^^, corresponding to 232 x 10^ BB pairs, was 
recorded at the T(45) resonance (center-of-mass energy 
Vi= 10.58 GeV). 

Charged particles from the e+e" interactions are de- 
tected and their momenta are measured using five layers 
of double-sided silicon microstrip detectors and a 40-layer 
drift chamber, both operating in the 1.5-T magnetic field 
of a superconducting solenoid. Photons, electrons, and 
hadronic showers from K'^ interactions are identified with 
a CsI(Tl) electromagnetic calorimeter (EMC). Further 
charged particle identification is provided by the specific 
ionization (dE/dx) in the tracking devices and by an in- 
ternally reflecting ring imaging Cherenkov detector cov- 
ering the central region. The steel of the magnetic-fiux 
return for the superconducting solenoid is instrumented 
with resistive plate chambers (instrumented flux return, 
or IFR), which are used to identify muons and hadronic 
showers from interactions. 

III. ANALYSIS METHOD 

A. Daughter candidate selection 

We reconstruct K'^ candidates through the decay 
TT+TT" only. We begin by forming all oppositely- 
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charged combinations of reconstructed tracks in the 
event. The invariant mass is required to be in the range 
of 490 to 506 MeV/c^, assuming the tracks arc pions. In 
addition, the probabihty of the vertex fit must 
be greater than 1% and the transverse component of the 
decay length is required to be greater than 8 mm. Fi- 
nally, the angle between the flight direction and the 
momentum vector must be less than 0.2 radians. 

We reconstruct candidates by identifying clusters 
of energy in the EMC and hits in the IFR that are iso- 
lated from all charged tracks in the event. The can- 
didates based on IFR clusters must have hits in at least 
two of the 16 to 19 layers in the detector. The can- 
didates detected by a hadronic shower in the EMC must 
have a calorimeter energy of at least 200 MeV, where this 
energy is from interpreting the calorimeter signal as an 
electromagnetic shower. Clusters in the EMC that are 
consistent with photons from 7r° decay are vetoed by 
two methods. If the candidate EMC cluster forms 
an invariant mass in the range of 100 to 150 McV/c^ 
with another EMC cluster with a calorimeter energy of 
at least 100 MeV, it is rejected. If the candidate clus- 
ter has a calorimeter energy greater than 1.0 GeV and 
two local maxima (or two "bumps" ) in the spatial distri- 
bution of the energy within the cluster, it may be from a 
high-energy tt*^, where the electromagnetic showers from 
the two photons arc merged into one cluster. If the two 
bumps within the cluster form an invariant mass greater 
than 110 MeV/ c^, the candidate is rejected. 

Further background rejection for EMC candidates 
is achieved by using a neural network trained on signal, 
qq continuum (where q = u,d, s, c), and BB Monte Carlo 
samples to distinguish K'^ clusters from fake clusters. 
The neural network inputs are the EMC cluster energy 
and the following six shower-shape variables; 

• The lateral moment LAT of the shower energy de- 
position [14] defined as LAT = Y.7=3 Eirf/{Eirl + 
-^27'o + X]r=3 -^i^i ) where the n crystals in the EMC 
cluster are ranked in order of deposited energy {Ei) 
, ro = 5 cm is the average distance between crys- 
tal centers, and is the radial distance of crystal 
i from the cluster center. 

• The second radial moment of the shower energy 
deposition, defined as ^ ■ Eirf where ri is the radial 
distance of crystal i from the cluster center. 

• The energy sum of a 3 x 3 block of crystals, centered 

on the single crystal with the most energy, divided 
by the larger 5x5 block, also centered in the same 
way. 

• The energy of the most energetic crystal in the clus- 
ter divided by the energy sum of the 3x3 crystal 
block with the most energetic crystal in the center. 

• The Zernike moments [15] ^2,0 and ^4^2 defined 
below. 



The Zernike moment An,m is defined as 



with 



1 



'tot 



and 



/4 2(a;) =^x^ -?>x^ 



where and <l)i are the radial and angular separation of 
crystal i with respect to the cluster center, E^ot is the 
total cluster energy, and i?o is a cutoff radius of 15 cm. 
Only candidates that pass an optimized cut on the 
neural network output are retained. This cut has a signal 
efficiency of 85% and rejects 70% of the EMC fake Kl 
background. 



B. B" candidate selection 

We reconstruct i?" candidates from selected if" clus- 
ters and pairs of K^^ candidates that do not share any 
tracks. We require the sum of the i^g momentum 
magnitudes in the center-of-mass frame to be at least 
2.1GeV/c, which ensures consistency with the kinemat- 
ics of a i?° ^ KKK decay. Only the direction of the 
is reconstructed, from the vector defined by the primary 
vertex and the center of the neutral cluster. We com- 
pute the momentum by constraining the K'^K^jK^ 
invariant mass to be the known mass of the B^ [16]. 
We reject fake candidates by using the difference be- 
tween the K1 transverse momentum, computed from the 
B^ mass constraint, and the transverse momentum along 
the K^l direction that is missing from the event. The re- 
constructed event missing transverse momcntiim is cal- 
culated without using the cluster in the momentum 
sum and projected along the computed direction of the 
K^. This missing transverse momentum difference (re- 
constructed minus calculated) is required to be greater 
than —0.5 GeV/c. This requirement and the previously 
mentioned requirement on the EMC neural network 
output were simultaneously optimized to give the great- 
est signal significance S/ \/S + B, where S and B are the 
expected number of signal and background events, as- 
suming a signal branching fraction of 5 x 10~^ [7]. 

The difference in energy AE between the reconstructed 
candidate and the beam energy in the center-of-mass 
frame is the main variable used to distinguish prop- 
erly reconstructed signal events from combinatoric back- 
ground. The missing transverse momentum difference 
distributions for the signal Monte Carlo sample and the 
background-dominated AE sideband {AE > 0.010 GeV 
in data) are shown in Figure 1. 

We distinguish the non-resonant three-body i?*^ decay 
from two-body decays to the same final state in this 
analysis. There are four two-body B° decays through 
charmonium that can go to the final state K'^K';!k1, 
namely: J/i^K^ XcoK, Xc2K°, and V(2^)^°- These 
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FIG. 1: Missing transverse momentum difference distribu- 
tions for the signal Monte Carlo sample (solid) and the 
background-dominated AE sideband {AE > 0.010 GeV) af- 
ter all other selection criteria (see text) have been applied. 
The signal Monte Carlo sample distribution has been arbi- 
trarily normalized to the number of entries in the sideband 
distribution. 



FIG. 2: Neural network output (NN) distributions for 
the signal Monte Carlo sample (solid) and the background- 
dominated AE sideband {AE > 0.010 GeV) after all other 
selection criteria (see text) have been applied. The signal 
Monte Carlo sample distribution has been arbitrarily normal- 
ized to the number of entries in the sideband distribution. 



are particularly unwanted, since they are from a color- 
suppressed tree decay amplitude, not the b sqq pen- 
guin decay amplitude. The Xc modes have unknown B'^ 
branching fractions. We veto candidates consistent 
with J5° XcoKI or B° Xc2K^^ where the Xc decays 
to K'^K^, by removing candidates with a K^K^ invariant 
mass in the range 3.400 to 3.429 GeV/c^ or 3.540 to 3.585 
GeV/c^ respectively. The combined _B" and daughter 
branching fractions for the J/tpK'^ and tp{2S)K'^ modes 
are 0.062 and 0.016 times 10~^ respectively. We expect 
about two events from two-body B*^ decays through char- 
monium in our final sample. 

We also remove B'^ candidates consistent with B'^ 
4>K1 with (f> K^K^ by requiring the invariant mass 
of both K^K^ combinations to be above 1.049 GeV/ c^, 
though as a cross check, we measure the branching frac- 
tion in only the (t)K^ region, where one K^K^ com- 
bination has an invariant mass less than 1.049 GeV/c^. 
These two-body B^ decay vetos are 96.6% efficient for 
the B^ K°K°K° signal Monte Carlo sample gener- 
ated with a uniform true Dalitz distribution. 



C. Event selection 

The main source of background is from random 
K'-^K'-^K'-^ combinations from continuum qq events. We 
combine 3 variables in a neural network that is trained 
using Monte Carlo samples to distinguish signal from 
continuum events. The first input is the cosine of the 
angle of the B^ momentum with respect to the beam 
axis in the center-of-mass frame 9g, which is flat for 
continuum background whereas the signal probability is 
proportional to sin^ 9*^ . The other two inputs are topo- 
logical variables that are commonly used to distinguish 
jet-like continuum events from the more isotropic parti- 
cle distributions in BB events. The first is the cosine 



of the angle 9t between the thrust axis of the can- 
didate in the center-of-mass frame and that of the rest 
of the charged tracks and neutral calorimeter clusters in 
the event. The distribution of | cos 0^1 is sharply peaked 
near 1.0 for combinations drawn from jet-like qq pairs, 
and nearly uniform for meson decays. The second 
is a linear combination of the zeroth and second angu- 
lar moments Lo,2 of the energy flow about the thrust 
axis. The moments are defined by Lj = X^iP* ^ |cos6'i|-' , 
where 6i is the angle with respect to the B° thrust axis 
of track or neutral cluster i, pi is its momentum, and the 
sum excludes the B^ candidate daughters. Distributions 
of the neural network output NN for the signal Monte 
Carlo sample and the background-dominated AE side- 
band are shown in Figure 2. We require NN > 0.5 to 
remove events that have little probability of being signal. 

After all of the requirements stated thus far, if an event 
has more than one B^ candidate, we choose the best one 
by selecting on the quality of the cluster. If there are 
one or more EMC if^ candidates, the best if ° candidate 
is the one with the highest cluster calorimeter energy. If 
there are no EMC candidates, the best candidate 
is the one with the highest number of IFR layers with hits 
in the cluster. If there is more than one K'^ pair that 
uses the same (best) A'" cluster, the best B^ candidate 
is the one with the lowest K'^ mass defined below 



X 



5rni 



+ 



Sm2 



where 6mi is the difference between the reconstructed 
invariant mass of K'^ candidate i and the known 
K'^ mass [16] and a is the invariant mass resolution 
(2.9 MeV/c2). 

Our final analysis sample contains 5892 events with 
AE in the range of -0.010 GeV to 0.080 GeV. Signal 
events appear mostly in the range of —0.010 GeV to 
0.010 GeV. The AE > 0.010 GeV region is dominated 
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by combinatoric background. The signal efficiency varies 
from 4% to 14% depending on the position in the Dahtz 
plot. The signal efficiency as a function of Dalitz plot 
position and a histogram of the reconstruction efficiency 
for 322 uniformly distributed points in the Dalitz plot 
are shown in Figures 3 and 4 respectively. For a uniform 
true Dalitz distribution, the average signal efficiency is 
8.1%. 




15 20 , . 

m2„^„(GeVV) 



in the e+e~ and ^jl^ n~ channels. We apply the same 
selection, projected missing transverse momentum differ- 
ence, and NN criteria as for our K2:K2;K^ selection. The 
ratio of the K^s^^Kl and J/ij^Kl Monte Carlo efficien- 
cies for these criteria is ^k°k%ki/^j/iPk1 — 0.96. We 
compare the number of ffited J/ipK^ events to the pre- 
dicted number of events, based on the known branching 
fractions and the Monte Carlo efficiency. The J/ ipK^ re- 
construction efficiency for all selection criteria is 12.5%. 
We find A^obs — 1420 ±56 J/^pK^ events, consistent with 
the predicted yield of A^o — 1426 ± 86. We correct the 
K'^K'^K^ efficiency by multiplying the MC efficiency by 
a correction factor of 0.96 ± 0.08, which is defined as 
(A^obs/A^o) X {^K%K%Kll^j/i,Kl)- The uncertainty on the 
correction factor includes the uncertainties on the rele- 
vant branching fractions for J/tpK'^, the statistical error 
from the J/ipKl fit, and our estimated uncertainty on 
relating the J/ipK^ and K^K'^K^ selection and recon- 
struction efficiencies. 

The AE resolution is better for EMC candidates be- 
cause the position of the hadronic shower is more pre- 
cisely measured. We also use the J/tpKl to measure the 
signal AE resolution separately for K'^ candidates recon- 
structed in the EMC and IFR. Figure 5 shows the fitted 
AE distributions for EMC and IFR K° candidates. 



FIG. 3: The reconstruction efficiency for signal as a function 
of position in the symmetrized Dahtz plot. 
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FIG. 4: A histogram of the reconstruction efficiency for signal 
for 322 uniformly distributed points in the Dalitz plot. 




0.08 



0.06 0.08 
AE,FR(GeV) 



D. J /ipKl control sample 

We use a sample of J/tpKl decays from the 

data to calibrate the reconstruction and selection effi- 
ciency and the AE resolution. The J/ip is reconstructed 



FIG. 5: Fitted AE distributions of the J/ipKl control sample 
for EMC (top) and IFR (bottom) Kl candidates. The points 
with error bars are histograms of the data sample. The solid 
curve is the total PDF. The dot-dashed curve is the signal 
PDF and the dashed curve is the background PDF. 
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E. Signal yield determination 

We use an extended unbinned maximum likelihood fit 
to determine the number of signal events in our final sam- 
ple. The likelihood for an event is the product of proba- 
bility density functions (PDFs) for the two main discrim- 
inating variables in the analysis: AE and NN. Separate 
PDFs are used for EMC and IFR candidates due to the 
difference in AE resolution and the fact that some back- 
ground channels mostly produce fake signal for EMC 
candidates only. The EMC and IFR samples are fit- 
ted simultaneously and the relative fraction of EMC sig- 
nal events is taken from the J/tjjK'^ control sample. Sep- 
arate PDFs are also used for signal, combinatoric back- 
ground, and three classes of background from B decays 
similar to our signal ("peaking backgrounds" described 
below in Section III F) . The combinatoric background in- 
cludes random ifg/fgif^ combinations from BB events 
as well as from continuum events. 

Figure 6 shows the AE and NN PDFs for all five 
fit components. The numbers of signal and combina- 
toric background events are determined from the fit. The 
peaking background yields are fixed to values based on 
known and estimated branching fractions and then varied 
in the evaluation of systematic uncertainties. 

The functional form of the signal PDF for AE is a 
triple Gaussian distribution. The mean and width of the 
core Gaussian distribution are determined separately for 
EMC and IFR X° candidates from the J/ipKl control 
sample and held fixed in the K^K^K^ fit. The remain- 
ing Gaussian parameters for the signal AE PDF are de- 
termined from the signal Monte Carlo sample and held 
fixed. The signal PDF for NN is a 4th-order polynomial. 
The polynomial coefficients are determined from the sig- 
nal Monte Carlo sample and held fixed in the fit. The 
combinatoric background AE PDF is an ARGUS func- 
tion [17]. The NN PDF for combinatoric background is 
the sum of a Ist-order polynomial and an ARGUS func- 
tion. The PDF shape parameters for the combinatoric 
background component are free in the fit. 



F. Backgrounds from other B decays 

Backgrounds from B decays to final states similar to 

K^K^KI can look similar to our signal in our discrimi- 
nating variables AE and NN . We call events from these 
decays "peaking backgrounds." The largest single source 
of peaking background is the decay B^ K^K^K'^. 
One of the Kg can decay to tt^tt^, where one or more 
of the photons from either tt" fakes the K'^ cluster in the 
EMC. The world average K°K°K° branching fraction 
from the Heavy Flavor Averaging Group (HFAG) [6] , [5] 
is B{B° ^ K^K^.K';^) = (6.2 ±0.9) x 10"^ We expect 20 
KgK^K^ events in our fit sample and vary this number 
by ±5 events in the evaluation of the systematic errors. 

Decays of the type KKK* can produce a /fgA'^/f^Tr 
final state, which can look similar to our signal if the 



momentum of the additional tt is low. The branching 
fractions for the relevant KKK* have not been mea- 
sured. We assume they are each of the same order as 
our expected signal, based on comparing the relative 
branching ratios of S — » Kit decays with B K*Tr and 
B Kp decays. We estimate a total branching fraction 
of 30 X 10~^ for this inclusive final state and vary this 
assumption by ±15 x 10~^ in the evaluation of system- 
atic errors. This gives us an expected 70 ±35 events from 
KKK* in the fit sample. In the evaluation of the 90% 
C. L. upper limit on the branching fraction, we set the 
KKK* event yield to zero. This is the most conservative 
assumption, since the signal and KKK* event yields are 
anti-correlated in the fit. 

Finally, the decay B^ K^KIti^ can also mimic our 
signal if the tt" is misidentified as a K^ in the EMC. 
The branching fraction for this mode has also not been 
measured. We estimate the branching fraction, relative 
to our signal, by assuming that the tree-to-penguin ratio 
for the three-body decays is the same as for B^ tt^tt" 
vs B^ Ar°7r° and that our signal branching fraction 
is about 5 x 10~^ [7]. Our estimate for the branching 
fraction is B{B" K'^^K'^^tt^) w 1.6 x 10"^ and we vary 
this assumed branching fraction in the range [0.1,5.0] x 
10~^ in the evaluation of systematic errors. This gives us 
an expected 2t_2 events from K^K^n*^ in the fit sample. 

The AE and NN PDFs for the three peaking back- 
ground components above are determined from Monte 
Carlo samples. In general, the AE shape peaks near 
AE = 0, though there is a substantial tail that extends 
to high AE values. We use the sum of an ARGUS func- 
tion and a Landau function for the peaking background 
AE' PDF. The NN shape is quite similar to the signal 
shape for all three peaking backgrounds. The form of the 
NN PDF is a polynomial. The AE and NN PDFs for 
the three peaking background components are shown in 
Figure 6: c, d, and e. 



IV. RESULTS 

Table I lists the results of the maximum likelihood fit. 

We find 23^^3 ± 6 events for the K'^K'^^K^ signal yield, 
where the first errors are statistical and the second error 
is systematic. The maximum likelihood fit bias of +0.3 
signal events was evaluated from an ensemble of data sets 
composed of fully simulated signal and B background 
Monte Carlo events. The combinatoric background for 
these datasets was generated from the PDF parameters of 
the fit to the data. This average bias of 0.3 events and the 
expc;cted 2.1 events from decays through charmonium 
are subtracted from the fitted signal yield. 

The systematic errors on the fitted signal yield and 
the signal branching fraction are listed in Table II. The 
additive contributions to the systematic error come from 
the uncertainty on our estimation of the KKK* com- 
ponent normalization (5.2 events), varying the yields of 
the other fixed peaking background components and the 
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FIG. 6: Distributions of AS (left) and iVA^ (right) for the 
five components of the fit: (a) K%K%K^ signal, (b) com- 
binatoric background, (c) KIkI-k° , (d) K^KIkI, and (e) 
KKK* . The combinatoric NN distribution is from the data 
sideband (Ai? > 0.010 GeV). The rest of the distributions are 
from Monte Carlo samples. 



fixed PDF parameters (2.6 events), our uncertainty on 
the fit bias correction (0.2 events), and the uncertainty 
on the charmonium background subtraction (1.1 events). 
The multiplicative contributions are from the uncertain- 
ties on the and reconstruction efficiencies (6.8 and 
8.0%), the number of BB events in the sample (1.1%), 
and the K^g —>■ tt+tt^ branching fraction. 

Figure 7 shows the fitted distributions of AE and NN, 
where a cut has been made on the fit variable not shown 
(e.g. there is a NN cut applied for the AE plot) to 
enhance the signal. The signal efficiency of this cut is 57% 



for the AE distribution and 71% for the NN distribution. 
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FIG. 7: Distributions of AE (top) and iVA^ (bottom) for the 
KsKsKI fit. The plot of AE (NN) is for events passing a cut 
on A'^A'" (AE) which enhances the signal. The signal efficiency 
of this cut is 57% for the AE distribution and 71% for the 
A' TV distribution. The points with error bars are histograms 
of data samples. The solid curves are total PDFs. The dashed 
curves are combinatoric backgrounds. The dot-dashed curves 
are peaking B backgrounds. The dotted curves are signal 
PDFs. 



The true distribution of K'^K^K^ events in the Dalitz 
plot is unknown. We use the average signal efficiency 
assuming a uniform true Dalitz distribution (8.1 %) to 
compute the ifgXgX^ branching fraction. We found 
no significant dependence of the signal AE resolution or 
the signal NN shape on the Dalitz plot variables. With 
these assumptions, we find a K^K^K^ branching frac- 
tion, excluding the 4> resonance, of B{B'^ K^K'^KD — 
(2.4^2 5 ± 0.6) X 10~^, where the first error is statisti- 
cal and the second error is systematic. The dominant 
systematic error is from the uncertainties of peaking B 
backgrounds (23% relative). 

Figure 8 shows a scan of the negative log likelihood as 
a function of the K'^K'^KI branching fraction, where the 
minimum negative log likelihood has been subtracted. 
In order to remove any dependence on our estimate of 
the total KKK* branching fraction, we conservatively 
fix the KKK* yield to zero in the scan of the log likeli- 
hood. We compute a one-sided Bayesian 90% confidence- 
level upper limit on the branching fraction assuming a 
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uniform prior probability for positive (physical) branch- 
ing fraction values. Systematic errors are included by 
convolving the fit likelihood with a Gaussian distribu- 
tion with a width corresponding to the total systematic 
error, excluding the uncertainty on the KKK* yield, 
since it is fixed to zero in the likelihood scan. The re- 
sult for the non-resonant three-body branching fraction 
is K°K^K^) < 7 Ax 10^^ at 90% C.L. Assuming 

the worst-case true Dalitz distribution, where the signal 
is entirely in the region of lowest efficiency (4%), the 90% 
C. L. upper limit on the branching fraction is 16 x 10^^. 




FIG. 8: Scan of the fit negative log likelihood as a function of 
the K^KsKI branching fraction, where the minimum nega- 
tive log likelihood has been subtracted. A uniform true Dalitz 
distribution has been assumed for the K'^KgK'l signal. The 
solid curve includes the systematic uncertainty with yield of 
KKK* fixed at to calcuate upper limit. 



TABLE I: Results of the fit for yields, branching fraction cal- 
culation, and branching fraction upper limit (UL) calculation. 
The K'^KlK'l efliciency below assumes a uniform true Dalitz 
distribution. The maximum likelihood fit bias and charmo- 
nium background are subtracted from the signal yield branch- 
ing fraction calculation. 



Mode 


K°K°Kl 


SK^l 


Events to fit 


5892 


210 


Signal yield 


23"'' no ± 6 


8 3"''^^ 

"•^—4.5 


Comb. Bkg 


5777 ± 79 


202 ± 15 


KKK* Bkg 


70 (fixed) 


NA 


KaK^Ki Bke 


20 (fixed) 


NA 


K'iK'^Tv^ Bkg 


2.3 (fixed) 


NA 


Estimated fit bias (evt) 


-1-0.3 


0.0 


ccK^s /I, 


2.1 


NA 


MC e (%) 


8.7 


6.1 


K° corr. (%) 


96.2 




Kl corr. (%) 


96.1 




Corr. e (%) 


8.1 


5.6 


ne, (%) 


47.5 


16.2 


Corr. e X n (%) 


3.8 


0.91 




2.4^11 ±0.6 




Stat, signf (a) 


1.0 


2.2 


Signf. w/ syst. (a) 


0.9 




90% CL UL B (xlQ-^) (stat.) 


6.3 


NA 


90% CL UL B (xlO-*) (incl. syst.) 


7.4 


NA 



As a cross check on the analysis, we have performed 
the fit only in the (f>K'^ region of the Dalitz plot where 
the invariant mass of one of the K'^K'^ pairs is less 
than 1.049 GeV/c^. The results are given in the sec- 
ond column of Table I. We find S.S^I^s signal events, 
which corresponds to a branching fraction of B{B^ — > 
(jjK^) = (4.0^2 2) ^ 10^®, where the errors are statistical 
only. This is consistent with the world average value of 
0.5 • BiB° ^ <j)K") = (4.3t^-]) X 10-6[16]. 

We checked our estimation of the KKK* peaking 
background yield by allowing it to float in the fit. This 
fit gave a KKK* yield of —54 ± 170 events, which is 
consistent with our estimation of 70 ± 35 events. 



V. SUMMARY 

We have searched for the decay B'^ K'^K'^KI using 
232 million BB pairs recorded by the BABAR experiment. 
We find no significant evidence for this decay. The central 
value for the branching fraction, assuming a uniform true 
Dalitz distribution for the signal and excluding the (j) 
resonance, is B{B° K^^K^K^) = (2.4j:|^ ± 0.6) x 
10~^, where the first error is statistical and the second 



error is systematic. This corresponds to a Bayesian 90% 
C.L. upper limit of B{B° K^^K^Kl) < 7.4 x lO'^. 
Assuming the worst-case true Dalitz distribution, where 
the signal is entirely in the region of lowest efficiency, the 
upper limit on the branching fraction is 16 x 10~^. Our 
results show that the K'^K'^K'l channel will be of limited 
use in understanding the b — > sqq penguin CP anomaly, 
due to the low efficiency times branching fraction, which 
limits the yield of signal events. 
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TABLE II: Estimates of systematic errors. Multiplicative sys- 
tematic errors are in percent while additive systematic errors 
are in events. The fit yield systematic error is due to fixed 
fitting parameters, mainly from K'^K'^KI AE PDF param- 
eters, and the uncertainties of fixed peaking B background 
yields. The fit bias error is one half of the bias. The ccKs/l 
error is the uncertainty of charmonium background subtrac- 
tion, estimated as one half of the subtraction. 



Quantity K°K°K° 

Additive errors (events) 

KKK* normalization 5.2 

Fit yield 2.6 

Fit bias 0.2 

ccKs/L 1-1 

Total additive (events) 5.9 



Physics (China), the Commissariat a I'Encrgie Atom- 
iquc and Institut National do Physique Nucleaire et de 
Physique des Particules (France) , the Bundesministerium 
fiir Bildung und Forschung and Deutsche Forschungsge- 
meinschaft (Germany), the Istituto Nazionale di Fisica 
Nucleare (Italy), the Foundation for Fundamental Re- 
search on Matter (The Netherlands), the Research Coun- 
cil of Norway, the Ministry of Science and Technology of 
the Russian Federation, and the Particle Physics and As- 
tronomy Research Council (United Kingdom). Individu- 
als have received support from CONACyT (Mexico), the 
Marie-Curie Intra European Fellowship program (Euro- 
pean Union), the A. P. Sloan Foundation, the Research 
Corporation, and the Alexander von Humboldt Founda- 
tion. 



Multiplicative errors (%) 

Kg efficiency 6.8 

Kl cfliciency 8.0 

Number BB 1.1 

BiB° ^ Kl n+n-) 0.4 

Total multiplicative (%) 10.6 

Total errors [^(10"'^)] 0.6 
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